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Abstract 
Using the variational method, a detailed study of the lowest m — 0,-1 electronic states of the 
exotic molecular ion H^'^ in a strong magnetic field, in the linear symmetric configuration parallel 
to the direction of the magnetic field is carried out. A extended study of the lag ground state (J.C. 
Lopez and A. Turbiner, Phys. Rev. A 62, 022510, 2000) was performed obtaining that the potential 
energy curve displays a sufficiently deep minimum for finite internuclear distances, indicating the 
possible existence of the molecular ion H^^ , for magnetic fields of strength B > 3 x 10^"^ G. It 
is demonstrated that the excited state l7r„ can exist for a magnetic field B — 4.414 x 10^'' G 
corresponding to the limit of applicability of the non-relativistic theory. 

Keywords: Strong Magnetic Field; Molecular Ion Hf^; Excited States. 






O 

Oh, 



oo ■ 

o 

O' 

O 

o 

Oh 



X 



PACS numbers: 31.15.Pf,31.10.+z,97.10.Ld 



I. INTRODUCTION 

The neutron stars were predicted by L. Landau jlj in 
1932. This kind of stars are characterized by an enor- 
mous magnetic field {B ~ 10^^ G -lO^^G). Since the 
discovery of neutron stars in 1967 the behavior of the 
molecular and atomic systems in such strong magnetic 
field has been object of intense investigations. The first 
studies were made in plasma physics and astrophysics by 
Kadomstev-Kudrayavstev [J and Ruderman Q , respec- 
tively. The results demonstrated that the binding energy 
of a molecular system grows substantially with magnetic 
field increase, while the size of the systems decreases. 
As a consequence of such a behavior, the formation of 
a molecular chain aligned in the magnetic field direction 
arise. 

Recent studies of the Coulomb systems in a strong 
magnetic field have been encouraged by the observation 
of the neutron star 1E1207. 4-5209 {Chandra X-ray obser- 
vatory 2002). These observations displayed two absorp- 
tion features at 0.7keV and 1.4keV. The origin of such 
lines could be associate to the presence of exotic molec- 
ular systems in the atmosphere of the above mentioned 
star 4] . In particular, the studies of one-electron molecu- 
lar systems, traditional and exotic have been intensified, 
such as H+ Hl+, Hl+ , {HeHf+, Hel+ y Lil+ (see for 
example p)- 

The focus of this article is a quantitative analysis of 
the lowest electronic states of the molecular system ij|"'" 
placed in a strong magnetic field. The variational method 
in the framework of non-relativistic theory is used, i.e. 
B < 4.414 X 10^^ G (Schwinger hmit). Specifically, an 
extension of the preliminary study p] for the ground 
state is made, and the possible existence of excited states 
is carried out. 

Atomic units are used throughout [h = e = rUg = 1) 
and the energy is expressed in Rydbergs (Ry). Spectro- 
scopic convention is used to denote the states. A num- 
ber indicates the principal quantum number, Greek let- 



ter (T, TT, (5 . . . indicate the magnetic quantum number for 
mi = 0, 1, 2 . . . respectively. The spatial parity is defined 
by g (gerade) for states with even parity and u (unger- 
ade) for odd parity states. Particularly, the ground state 
is indicated as Ictq. 



II. THE MOLECULAR ION Hl+ 



The exotic molecular ion Hf^ is a system constituted 
by four protons and one electron. For the present work 
the protons arc consider over a straight line (linear config- 
uration) parallel to the lines of magnetic field B (parallel 
configuration) . The magnetic field B is taken as constant 
and homogeneous in z direction. The Born-Oppenheimer 
approximation of zero order is assumed (protons are con- 
sider infinitely massive). The positions of the four pro- 
tons over z axis are denoted as Ri, R2, — R3 and — R4. 
We will suppose also that Ri = R3 and R2 — R4 (sym- 
metric configuration (see Fig^). 

In order to describe the magnetic field B = (0,0,-6), 
the vector potential in the symmetric gauge A — 
-|-(— y, a;, 0) is assumed. The Hamiltonian that describes 



the molecular ion H^ (in cylindrical coordinates) is 



H 



LB 



B^p^ 



+ V{p, z) , 



(1) 



y^ is the distance from the electron to 



where p — x 

the z axis. 4 is the z component of the angular mo- 
ment (which is a integral of motion with eigenvalue m 
-magnetic quantum number). V(p, z) is the Coulomb po- 
tential of the interaction between the electron and each 
charged center, plus the classic repulsion term between 




-R3 
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FIG. 1: The molecular ion H'^'^ in a linear configuration 
placed in a constant and homogeneous magnetic field directed 
according to z axis, B = (0,0,-B). The protons positions are 
indicated by the four points over z axis. The position of the 
electron is marked by a point out of the axis. 



the charged centers. The potential is given by 
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where r^ = ^J p^ + (2; — Zi)"^ is the distances between the 
electron and the corresponding proton {zi are the coor- 
dinates of each proton over the z axis) . 

In the Hamiltonian (^ the spin contribution is ne- 
glected because it only changes a reference level for the 
total energy. 



III. THE TRIAL FUNCTIONS 

In order to study the system H\'^ the variational 
method is used. For the choosing of the trial functions, 
the fitting criteria used by Turbiner-Lopez was followed 
(see for example @). These criteria have proved to give 
much precise results, qualitatively and quantitatively, for 
different Coulomb systems in strong magnetic fields. Ba- 
sically these are: 



1.- 



_ Mt 



for which the trial function 

-A + T4,with 



The potential Vt — ^ 

il^t is an exact eigenfunction of Tit 
zero energy, should reproduce the original potential 
near the Coulomb singularities and the harmonic 
oscillator behavior at large distances. 

2.- The symmetries of the system must be included. 

3.- For the ground state, the trial function should not 
vanish inside the domain where the problem is de- 
fined(Theorem of Perron). 



In addition, since the charged centers (protons) are 
identical, the trial function is chosen symmetric with re- 
spect to the interchange of these, i.e., invariant under 



transformations P^ (r^ 



/), «,J = !■ 



The trial function is a linear combination of terms that 
fulfill the requirements above mentioned. The general 
form of the ansatz is 



(3) 







*m = ZYm ^0 , 
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e-'^«T,(4a) 



(4b) 



Here ^^ corresponds to the general trial function for the 
ground state (to = 0). ai.2,3.,4 and /3 are variational 
parameters. The sum is over all permutations of the 
parameters ai_2.3,4- Excitations for different values of 
the magnetic quantum number m are included with the 
Um term. The trial functions '^m given by Eq. (O cor- 
responding to different values of m are orthogonal and 
result suitable to describe the lowest states inside each 
family of own states with quantum number m defined. 



A. Gauge Rotation 

Considering the trial function ^„i given by Eq. (0) and 
the Hamiltonian H (see Eq. (^), the total energy (for a 
state with a quantum magnetic number to) is 



yp2|™l*A„*„rfr 



{as, (J si f^2\rn\^^. 



where 

hm =U-,^iiUm = Pm +Bm+ — p2 + V{p, z) , (6) 

is the gauge rotated Hamiltonian, and p^ = ^m^P ^m is 
the covariant moment . From here on, Et^^^ will denote 
simply as Et. This formulation is very suitable since 
each Hamiltonian h„i (for m constant) describes a family 
of eigenstates with the quantum number m. 

Moreover, the ground state and the excited states are 
described by trial functions which have the same struc- 
ture (given by the function H4al) '). The only difference 
appears as a weight factor p^'™'. 

The potential associated to an arbitrary trial function 
4>g is now determined by 



Vr 



(m) 









(7) 



where the Laplacian covariant A is defined as 

A = A + 2|m| 



ld_ 
pdp 



j_d_ 



(8) 



Hence, the general form of the potential associated with 
each trial function term '^g, (Eq. (|4a|l ) is llj 



V, 



(m) 



i—1 i,k—l 



(9) 



where 



Hi Hfe 



nrk 



-{P^ + {Z - Zi){z ~ Zk)} ; 



and hi = (r — ri)/|r — r^j is a unitary vector oriented 
according to the line which connects the ith position 
charged center, with the position of electron r. The po- 
tential exactly reproduces the original potentials of 
Coulomb for ai — a2 — a^ = a4 = 1, and of harmonic 
oscillator for /3 = 1. 



IV. TRIAL FUNCTION TERMS *p 

The election of the most general trial function ^^ 
(Eq. (|4a|) ') was in no way arbitrary. The a parameters 
may be interpreted as the effective charge of the protons 
and /3 parameter as the effective coupling of the electron 
with magnetic field B. With this characterization of the 
variational parameters, it is possible obtain different trial 
functions each one with certain physical interpretation. 

Let us consider for example the case in which the in- 
teraction of the electron with each one of the protons is 
coherent. That is, the electron interacts in the same way 
with each proton. To describe such physical situation, 
it is enough to take ai — a2 — a^ ~ a^ = ai in the 
Eq. H4a() . Defining this function as -0^^^ the following is 
obtained 



V'l 



1) 



g-Q!i(ri+r2+r3+r4) g-/3l 



This function is the Heitler-London type function (used 
to describe the system H2 without magnetic field) mul- 
tiplied by the lowest Landau orbital. The results of the 
calculus verify (a posteriori) that this function gives the 
major contribution in the description of the system near 
its equilibrium position. 

In order to describe the physical situation at large in- 
ternuclear distances, all a' s parameters are equal to zero, 
except one, which is denote as a2. The trial function is 



^r^ = {^ 



which corresponds to a Hund-MuUiken type function 
(used to describe the system H^ without magnetic field) 
multiplied by the lowest Landau orbital. This function 
describes the interaction of the electron with each of the 
protons independently {incoherent interaction). 

It is now clear that different degenerations from the 
general function (Eq. H4a|l 'l lead to different physical de- 
scriptions of the system, and therefore to several trial 
functions. In this case and in order not to discriminate 
between one and another of the different descriptions, the 
linear combination of all the possible particular cases was 
considered as trial function. Explicitly, the trial function 
used in the current study is 

All the terms included in the last equation admit a cer- 
tain physical interpretation. This are summarized in Ta- 
ble m where we also present the variational parameters 
number contained in each one of them. A detailed anal- 
ysis of each term appears in 7] . 



Function 
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5 


"H + H + H + H" 



TABLE I: The trial function terms *t (Eq.Hnj along with the 
number of variational parameters (N.V.P.) and its physical 
interpretation (see text). 

Another degeneration of the Eq. H4a|l happens when 
a's parameters are equal by pairs, ai = a2 — a; a^ — 
ai = a' , {a ^ a')- The physical picture correspond to 
an electron shared by two protons of different form. The 
situation may be interpreted as an interaction H^ + H^ 
between two Hydrogen molecular ions, even though of 
course we only have one electron (these type of physical 
interpretations are denoted between quotation marks on 
Table P). In the same way, in the most general case (all 
a's parameters are different), the corresponding function 
is a sum of products of four Hydrogen functions by one 
Landau function (see Eq. H4a|l 'l. This situation describes 



a mixed state of four Hydrogen atoms sharing only one 
electron. 

Along with the v4i_...^ii coefficients of the linear com- 
bination (Eq. Hl()|l ^ and the positions i?i,....4 (see Fig. QJ 
there is a total of 48 variational parameters. The normal- 
ization condition of the trial function allows to keep one 
of the A coefficients as constant. The symmetry restric- 
tion (i?3 — i?i,i?4 — R2) fixes two parameters. Finally 
the number of effective variational parameters is reduced 
to 45. 



Each step in a particular minimization consists in keep- 
ing some parameters as constant while others were varied. 
The election of which parameters must be kept constant 
and which free was always based on the physical inter- 
pretation of such parameters. In no way the calculation 
is automatic, rather, a sophisticated artisan treatment is 
required. With this procedure, a variational parameter 
configuration was obtained in which the systems total 
energy value is minimum. Each minimization procedure 
for a given magnetic field as a general statement took 
various months of calculations (see Q)- 



V. COMPUTER TECHNIQUE 

All calculations were made in a DELL PC with two 
Xeon processors of 2.8 GHz each. For minimization 
calculations the minimization routine MINUIT of the 
CERN-LIB library was used. Numerical integration is 
carried out using the routine DOIFCF from the NAG-LIB 
library. The relative precision in the numerical integra- 
tion was of 10""'^^. It is important to mention that once 
all variational parameters are determined, the calculation 
of the variational energy takes a few seconds. 



A. Minimization Strategy 



VI. RESULTS 

In this section the total energy Et, the binding en- 
ergy i^ Eb {Eh = B — Et), and the equilibrium distance 
of protons (i?3 = Ri, R4 ~ R2) for the ground state lag 
and the excited state l7r„ of the molecular ion H^^ are 
presented. The trial function used is shown in Eq. I|1U|I . 
For a comparison with the results obtained by Turbiner- 
Lopez 13j the calculations were carried out considering 
two different conversion factors to atomic units. These 
factors are Bq = 2.35 x 10^ G (the one adopted in this 
work) and Bq = 2.3505 x 10^ G ( the one adopted in ^). 



The variational energy calculation includes two techni- 
cal aspects: the numerical evaluation of integrals in two 
dimensions and minimization of energy with respect to 
the variational parameters. The objective in the mini- 
mization process is to find a point in the multidimen- 
sional parametric space that minimizes the total energy 
of the system (Eq. (O). Since we have 45 variational 
parameters, the search of the minimum is a formidable, 
complicated and tedious task. For this reason it is not a 
good idea to start the procedure making a minimization 
in which all parameters vary simultaneously. The strat- 
egy used consisted in making minimization sequences in 
the following way 

• Ansatz by Ansatz, knowing clearly that between 
the eleven terms that appear in the Eq. (|10|l . the 

most general V'po should give us the best (lowest) 
value in total energy. 

• Combinations of two by two Asantze, taking ipP^ as 
the common Ansatz to all these combinations. 

• Some combinations of three by three Ansatze, tak- 
ing as ground the previous combinations, where one 
more Ansatz was included. 

• Beginning from that combination of three Ansatze 
with the best total energy value, new Ansatze were 
incorporated until there was a total of eleven of 
them. 



A. The Ground State lag 

For the ground state lag each Ansatz presents a mini- 
mum in the energy curve for the finite parameter values 
Ri,...,4- ^^ Table ITT) the results are present. The results 
carried out by Turbiner-Lopez ij| (also present in Table 
ITT|) only consider the four terms that appear in Table Q] 

A total energy curve of the system as function of R2 
was obtained (see Fig. O. For this the parameter i?i 
is considered constant to the corresponding value of the 
minimum total energy configuration (in symmetric con- 
figuration). The i?3 and R4 parameters are now con- 
sidered as variational parameters. In other words the 
symmetry restriction is eliminated. The curve presents a 
well defined minimum for finite internuclear distances in 
the symmetric configuration. The height of the barrier is 
AE = Ejnax — Emin = 0.242 Ry. The asymptotic behav- 
ior of this curve for large i?2 shows that Ri+R^ w 0.1113 
(see Fig.^ and R4 — R3 « 0.1110. These values are very 
similar to the equilibrium distances of the ion H^~^ (see 
Table |IV||. This indicates that a decay hI~^ -^ Hf^ +p 
is possible. The energy as function of i?i , R2, R3 and R4, 
generates a surface. In particular the curve presented in 
Fig. 121 corresponds to the cut throughout the valley of 
this surface. 

The electronic distribution corresponding to ground 
state lag is shown in Fig. 13 We can observed that the 
electronic cloud is mainly concentrated around the inter- 
mediate point to the protons. 



B(G) 



5o(xlO'^G) ETJRy) EtjRy) Rija.u.) _R2(a.u.) 



3x 10^ 



2.35 



12727.5422 38.4153 0.0635 0.2763 Present 



2.35 
4.414 X 10^^ 2.3505 
2.3505 



18739.2526 43.7261 0.0574 0.2336 
18735.2602 43.7230 0.0574 0.2336 
18735.4675 43.5157 0.0572 0.2326 



Present 
Present 
Turbiner-Lopez " 



"Not published. 



TABLE II: Results obtained for ground state lag of the molecular ion H^ in the presence of a magnetic field B for different 
values of the conversion factor Bo and in comparison with those obtained by Turbiner-Lopez. 



Et "="' 

18740 



18737 
18736 
18735 
18734 
18733 
18732 
18731 




FIG. 2: Total energy curve of the molecular ion H^'^ in 
the ground state lag as a function of parameter R2, which 
presents a minimum for finite equilibrium distances. This 
in linear configuration and parallel to the magnetic field 
B = 4AU X 10^^ G. 



FIG. 3: Electronic distribution of the molecular ion lag in 
its ground state in a magnetic field B = 4.414 x 10^'^ G, with 
Bo ~ 2.35 x 10^ G, in linear and parallel configuration to the 
direction of the magnetic field, considering that Ri — R3 and 
R2 = R4 (see Fig. 1). 



B. The Excited State l7r„ 

The calculations were made for the two values of the 
conversion factor Bq and it is shown in Table UTTl Also 
for this state, each Ansatz presents a minimum in the 
energy curve for the finite parameter values i?i,...,4. 

For the excited state ItTu the electronic distribution is 
shown in Fig. 0] The positions of the protons are indi- 
cated by points in the plane (x, z). The fact that appar- 
ently there are only three points present and not four, it 
is due to the fact that the distances between the inter- 
nal protons is very small in comparison with the whole 
dimension of the system (see Table llllll . This electronic 
distribution presents seven maximums (Fig.^Jb)) located 
in neighboring regions to the protons. From these, those 
that surround internal protons (see Fig. ^ are ~ 10^ 
times more pronounced that the minors of them. A 
comparison of the total energy and the equilibrium dis- 
tances [iJI between molecular ions hI~^ and H2 , both 
for the excited state l7r„, shows that both are very simi- 
lar (see Table IrVJ) . indicating that if|+ -^ H^ + 2p could 
be a possible channel of decay of the molecular system 
HI~^ . If a comparison of the total energy of the systems 



ff|+ and i/|+ is carry out (Table HVJI, we found that 
the molecular ion H^^ could have a possible channel of 
decay H^^ — + -^^3^ + P- Among the systems H, H^ , 
Hf^ and Hf^^ the first of them, iJ, results as the one 
with the greatest total energy in the rank of the mag- 
netic fields — 4.414 x 10^^ G, for which a possible decay 
7? -I- 3p is discarded. 



TrZ+ 
«4 



Finally, it is important to mention that if we relax the 
symmetry conditions (i?i,...,4 are arbitrary) in the mini- 
mization, the system returns to the symmetric configu- 
ration. This is fulfilled for the ground state {m = 0) and 
for the excited state (m = —1). 



VII. CONCLUSIONS 

Applying the variational method, a detailed study of 
the lag and ItTu electronic states of the exotic molecu- 
lar ion Hf^ in a strong magnetic field was carried out. 
The magnetic field B, was considered to be constant and 
homogeneous. 



B(G) Bo(xlO''G) ET(Ry) 


E6(Ry) 


Ri(a.u.) 


R2(a.u.) 


4.414 X 10" 2.35 18742.8185 
2.3505 18738.8254 


40.1603 
40.1578 


0.0651 
0.0651 


5.5534 
5.5537 



TABLE III; Results for the excited state Ittu of the molecular ion H\ 
Bo — 2.35 X 10^ G in linear, parallel and symmetric configuration. 



3+ 
4 ■ 



For a magnetic field B = 4.414 x 10^^ G, with 





State 


ET(Ry) 


E6(Ry) 


Req(a.U.) 


Ht 


lag 


18728.477 


54.5018 


0.1016 




IcTu 


18750.07 


32.912 


2.021 




ItTu 


18741.89 


41.09 


0.130 




iTTg 


18757.273 


25.7054 


2.237 




I5g 


18747.527 


35.407 


0.148 




1<5„ 


18761.18 


21.80 


2.230 




2Gg 


18781.576 


1.402 


3.120 


Hl+ 


lag 
lcr„ 


18727.7475 
Not exist 


55.2312 


0.110 




ItTu 


18742.7564 


40.2223 


0.145 




iTTg 


Not exist 








iSg 


18748.9067 


34.0720 


0.167 




ISu 


Not exist 






Hl+ 


lag 


18739.2526 


43.7261 


0.1149 




ItTu 


18742.8185 


40.1603 


0.1302 



TABLE IV: Comparison between molecular systems with one 
electron: _ff^ Ig, H^'^ ^ and H^'^ (present) in different 
states. All of them in linear, symmetric and parallel configu- 
ration in the presence of a magnetic field B = 4.414 x 10^'^ G 
(Bo = 2.35 X 10^ G). In H^+ the equihbrium distance Re, 
corresponds to the distance between to adjacent protons. For 
the system H^'^ this was considered as the separation between 
the two internal protons Reg = 2Ri (see Fig. 0. 



A. Ground State Itjg 

An improvement in the variational energy of ~ 0.2% 
was obtained in comparison with the first calculations 



presents in the literature (Turbiner-Lopez). The equilib- 
rium distances in both case are very similar. This in the 
symmetric configuration. An important remark is that 
the system H^^ can already exist for a magnetic field of 
intensity B == 3 x lO" G. 



B. Excited State l7r„ 

The existence of the excited state of the molecular ion 
in a magnetic field B = 4.414 x 10^'^ G is seem pos- 
sible. This state presents a peculiarity in its equilibrium 
distances (see Table 1111}) . i?2 turns out to be ~ 80 times 
Ri . The fact that the total energy of the molecular ions 
Hf^, H^~^ and H2 is very similar, indicates that the ion 



h!+ 



Hi^ is very unstable in this state. In this sense it is very 
unstable since any small disturbance could originate a 
decay. For which the existence of other excited states 
does not seem possible, at least for this magnetic field 
value. 
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